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ABSTRACT 

The GallCS (Galaxies in Cosmological Simulations) semi- analytical model of hierar- 
chical galaxy formation is used to investigate the effects of different galactic properties, 
including star formation rate (SFR) and outflows, on the shape of the mass metallic- 
ity relation and to predict the relation for galaxies at redshift z = 2.27 and z — 3.54. 
Our version of GallCS has the chemical evolution implemented in great detail and is 
less heavily reliant on approximations such as instantaneous recycling. We vary the 
model parameters controlling both the efficiency and redshift dependence of the SFR 
as well as the efficiency of supernova feedback. We find that the factors controlling the 
SFR influence the relation significantly at all redshifts and require a strong redshift 
dependence, proportional to 1 + z, in order to reproduce the observed relation at the 
low mass end. Indeed, at any redshift, the predicted relation flattens out at the high 
mass end resulting in a poorer agreement with observations in this regime. We also 
find that variation of the parameters associated with outflows has a minimal effect on 
the relation at high redshift but does serve to alter its shape in the more recent past. 
We thus conclude that the relation is one between SFR and mass and that outflows 
are only important in shaping the relation at late times. When the relation is stratified 
by SFR it is apparent that the predicted galaxies with increasing stellar masses have 
higher SFRs, supporting the view that galaxy downsizing is the origin of the relation. 
Attempting to reproduce the observed relation, we vary the parameters controlling 
the efficiency of star formation and its redshift d ependence and compar e the predicted 
relations with lErb et~aTI (|2006f ) at z = 2.27 and lMaiolino et al.l (|2008t ) at z = 3.54 in 
order to find the best-fitting parameters. We succeed in fitting the relation at z = 3.54 
reasonably well, however we fail at z — 2.27, our relation lying on average below the 
observed one at the one standard deviation level. We do, however, predict the observed 
evolution between z = 3.54 and z = 0. Finally, we discuss the reasons for the above 
failure and the flattening at high masses, with regards to both the comparability of 
our predictions with observations and the possible lack of underlying physics. Sev- 
eral of these problems are common to many semi-analytic/hybrid models and so we 
discuss possible improvements and set the stage for future work by considering how 
the predictions and physics in these models can be made more robust in light of our 
results. 

Key words: methods: N-body simulations - galaxies: evolution - galaxies: haloes - 
galaxies: formation 



1 INTRODUCTION 

The existence of a relation between luminosity (or mass) 
and metallicity in irregular and blue compact galaxies 



was first prop osed by Lequeux et al.l (Il979l) and later 
confir med by ISkillman et al. I l| 19891 ). iGarnett fc Shields! 
ll 19871) extende d the relation to spiral galaxies. Recently 
iTremonti et al.l |2004l ) examined the relation at z ~ 
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using ~53000 local star-forming galaxies in the SDSS 
and found that 12 + \og(0/H) increases steeply for 
10 8 ' 5 < M*/M© < 10 10 ' 5 but flattens at larger masses. 

Several explanations for such a rela t ion h ave been put 
forward. For instance, iTremonti et al.l l|2004h suggest that 
galaxies form similar fractions of stars independent of their 
mass. Less massive galaxies are simply less efficient at 
retaining gas due to their shallower potential wells and lose 
newly produced metals to g alactic outflows. Another possi- 
bility (|Maiolino et al.ll2008l ) is that more massive galaxies 
had higher specific star formation rates (SFR per unit 
stellar mass) or a higher star formation efficiency (SFR per 
unit gas mass) in the past so form a larger fraction of stars 
in the same Hubble time. In this case the higher metallicity 
is due to the conversion of a larger amount of primordial 
gas into C, N and O. This scenario is known as galaxy 
down-sizing fo r which there is a pl e thora of observational 



evidence (e.g Perez-Gonzalez et all |2008l )). An alternate 



explanation |Koppen et alj 2007 ) is a variation in the initial 



mass function (IMF) in different star forming environments. 
It is not known at present which process is responsible for 
the ubiquitously observed relation however it is likely that 
all three contribute to some extent. Interesting ly, a similar 
relation holds for the st ars of gas-poor galaxies (|Faberll 19731 ; 
iBrodie fc Huchralll99ll ). 

A lot of effort has been made to understand whether an 
underlying common physical origin for the above mentioned 
relation exists, as well as to determine the highe st redshift 
at wh ich the mass-metallicity relation holds. lErb et all 
(2006) have measured the relation at z a 2.27 using 
the [Nil] /Ha ratio in stacked spectra for a sample of 
87 rest-fr ame ultraviole t -select ed star-forming galaxies. 
Recently, iMaiolino et al.l |2008T ) have used deep near-IR 
spectroscopy of H/3 and [OIII]5007 shifted into the K band 
as well as [Oil] 3727 and [NeIII]2870 shifted into the H band 
to measure the relation for nine star fo r ming galaxies at 
z ~ 3.54. More recently, iMannucci et alj (|2009h have used 
near-IR spectroscopy of the optical lines [OII]A3727, H/3, 
and [OIII]A5007 for a sample of 10 Lyman- break galaxies 
(LGBs) at z ~ 3 to derive their SFR, metallicity, gas 
fraction and effective yield. Using optical, near- IR and 
Spitzer/IRAC photometry, they have measured the stellar 
mass of each galaxy in order to guarantee a robust estimate 
of the mass-metallicity relation. Finally, IMannucci et al.l 
(|2010T) introduce the SFR-mass-metallicity relation and 
show that it does not evolve with redshift up to z ~ 2. They 
argue that the apparent evolution of the mass-metallicity 
relation inferred by past works at redshifts below z ~ 2 
is only due to selection effects, where independent surveys 
sample different areas of the SFR-mass-metallicity "surface" 
at different redshifts (in particular those with a higher SFR 
at higher redshifts). A strong evolution still occurs for z > 2. 

Previous attempts at modelling hierarchical galaxy evolu- 
tion through semi-analytic models and JV-body simulations 
have shown that a mass-metallicit y relation is predicte d 
also in the hierarchical scenario jDe Rossi et al.1 l2007h . 
In one such simulation |Pe Rossi et al.1 "( 20071 ) predicted 
the relation over the redshift range < z < 3; however 
the predicted metallicities are always larger than the 



observed ones. The authors attribute this to the lack of 
supernova feedback in t he model and thus t he lack of 
outflows from the galaxies. |Pe Lucia et al. ( 2004) predicts a 
relation at z « for three models characterized by different 
feedback processes. It was found that all three models 
predict galaxies whose average metallicities lie within one 
standard dev iation of the median m ass-metallicity relation 
observed by Tremonti et all J2004I ), at variance with the 
results of iFinlator fc Dave] (I2008T ). who found that the 
predicted relation depends strongly on the galactic outflow 
model used. They also found that a lack of outflows lead 
to galaxies that we re too enriched in me tals, consistent 
with the results of De Rossi et all |2007). According to 
IFinlator fc Pavel l|2008h only a model where the wind was 
momentum drive n could reproduce the observed relation 
l|Erb et al.l 120061 ) . In this simulation the relation was 
reproduced reasonably well within experimental uncer- 
tainties and could match the observed relation in slope, 
amplitude and scatter. By incorporating mass dependent 
galactic winds in to the parallel tre e -SPH code GADGET-2 
l|Springelll2005l ). iKobavashi et al. | l|2007h predicted a rela- 
tion at z s a 2 t hat is consistent with the observations of 
lErb et al.1 l|2006h for massive galaxies only. Their predic- 
tions exhibit a significant scatter and their model suffers 
from the fact that it does not predict the termination of 
star formation in massive galaxies at late times. Finally, 
using A r -body/hydrodynamical simulations iMouchine et al.l 
(2008) were able to reproduce the observed relations over 
the range < z < 1.2 with minimal scatter although at 
higher redshifts the simulated relation predicts galaxie s 
with higher metallicities than are observed l|Erb et al ]|2006l ). 



As noted in IMaiolino et al.1 (|2008l ). where we address the 
reader for a more extended discussion of the single cases, 
many of the above mentioned attemp ts have failed to re - 
produce the observed relation at z ~ 3. ICalura et al.l (|2009l ) 
have recently been successful in predicting the relation 
over the range 0.07 < z < 3.54 using a chemical evolution 
model that predicts the relation separately for galaxies 
of different morphological type. Rather than collectively 
fitting the relation they have found that the relation at 
low redshifts is best fit by considering only spiral and 
irregular galaxies whilst at intermediate redshifts (z ~ 2) 
the relation is best fit by a mixture of proto-spirals and 
proto-ellipticals. At z — 3.54 they predict that the relation 
is best fit by proto-ellipticals alone. In this work we attempt 
to reproduce and explain the observed relation at z = 2.27 
and z = 3.54 by testing aspects of the galactic outflows (in 
common with many of the above models) as well as other 
methods including different IMFs and a redshift dependent 
SFR. In order to do this, we make use of an up-to-date 
version of GallCS in which a detailed treatment of the 
chemical evolutio n has been implem ented; details of which 
can be found in iPipino et al.1 l|2009l ). This model - which 
tracks the evolution of H, He, O and Fe - has the chemical 
evolution implemented in great detail. It does not rely on 
the instantaneous recycling approximation but instead uses 
a self-consistent prescription for both type la and Type II 
supernovae ejecta, which includes the effects of finite stellar 
lifetimes. We consider this aspect the main innovation of 
our approach. 
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We assume a flat (critical density) lambda cold dark 
matter cosmology (ACDM) with cosmological parameters 
taken from the W MAP three year results (WMAP3) 
jSpergel et alj l2007h . These are Q = 0.24, Ov = 0.76, 
h = 0.73 and £Ib = 0.06h where all the symbols have their 
usual meaning. 

In section [2] we describe the GallCS model and introduce 
the free parameters that it uses to control the star forma- 
tion efficiency, redshift dependence of the SFR, supernova 
feedback and the efficiency with which gas that was once 
ejected by the galaxy is re-accreted. In section [3] we show 
the data with which we compare our predictions and dis- 
cuss their uncertainties, mainly due to calibration issues. In 
section |4] we determine the effect of varying the free param- 
eters and IMF upon the predicted relation and compare our 
predictions with observations. In section [5] we then present 
the results of using the best-fitting parameters for the en- 
tire population of model galaxies and interpret the results 
physically. In section [6] we draw our conclusions. 



2 THE GALICS MODEL 

GallCS is a semi analytic hybrid model of hierarchial 
galaxy formation that combines the output of large A^-body 
cosmological simulations to track the evolution of baryonic 
matter in galaxies through their dark matter haloes. The 
evolution of galaxies is tracked using halo merging trees 
that follow the hierarchial evolution of small objects at 
early times that may or may not develop into larger 
ones through merg ing processes or accretion of matter 
ijHatton et alJ 12003). GallCS assigns a morphology to a 
galaxy instantly after a merger based on the ratio of the 
bulge to disc B-band luminosities. In outline, as hot gas 
cools and falls into the centre of its dark matter halo, it 
settles in a rotationally supported disc. The galaxies remain 
pure discs if their disc is globally stable and they do not 
undergo a merger with another galaxy. In the case where a 
significant merger occurs, we employ a recipe to distribute 
the stars and gas in the galaxy between three different 
components in the resulting, p ost-merger galaxy, t he disc, 
the bulge, and a star-burst fsee lHatton et al1l2003l ). In the 
case of a disc instability, we simply transfer the mass of the 
gas and stars necessary to make the disc stable to the burst 
component, and compute the properties of th e bulge /burst 
in a similar fashion to that described in lHatton et alj 
(2003). The star-burst scale is rb urs t = 0.1rb u i ge , so that 
the characteristic timescale for the star formation is shorter 
than in the bulge, hence leading to a faster consumption 
of the gas. The star formation rates are even higher than 
those in the discs, but have an instantaneous duration. The 
burst stellar population becomes part of the bulge stellar 
population when the stars have reached an age of 100 Myr. 
Since we will only be using the model predictions for high 
redshift it does not make sense to classify the galaxies 
by local (z = 0) standards and so here we consider all 
of the predicted galaxies no matter the GallCS assigned 
morphology. 

The simulation models the universe as a box of co- 
moving length of 150 Mpc. Like any numerical simulation, 



GallCS has a finite baryonic mass resolution. The minimum 
baryonic mass that we consider resolved in this simulation 
is 2xlO 9 M , which is a factor of ten lower than the fiducial 
GallCS value (Hatton et al. 2003). The baryonic gas in 
galaxies is initially primordial comprising of Hydrogen and 
Helium. The metal content increases as time passes due to 
the synthesis of these elements in stars during their lifetime 
and their subsequent release into the inter-stellar medium 
(ISM) upon the stars death. A detail ed description of the 
entire GallCS model may be found in lHatton et alj (I2003T ) 
and an updated version (as far as the implementation of the 
chemical evolution of O and Fe is concerned) of the model 
that w e use in this paper can be found in IPipino et al.l 
(|2009h . 

The main n o velty of the present version of GallCS 
l|Pipino et alj [2009) is the implementation of a self- 
consistent treatment of the chemical evolution with finite 
stellar lifetimes and both type la (SNIa) and type II (SNII) 
supernovae ejecta. In practice, we follow the chemical 
evolution of only four elements, namely H, He, O and Fe. 
This set of elements is good enough to characterize our 
simulated galaxies from the chemical evolution point of view 
as well as small enough in order to minimize computational 
resources. Also, the reader should remember that O is the 
major contributor to the total metallicity. 

We adopt the yields from llwamoto et al. I <|l999h and 
references therein for both SNIa and SNII. The reader 
should note that a change in the stellar yields will introduce 
a systematic offs et of a few tenths o f a dex in the model 
predi ctions (see iThomas et alj 120071 ; IPipino fe Matteuccil 
1 2004) . hence it might leave room for some fine-tuning for 
a suitable choice of stellar nucleosynthesis. However, being 
only an offset, this change cannot create nor modify the 
slope of the predicted relations. But, most importantly, 
the successful calibration of our m odel with element ratios 
observed in Milky Way stars (see IPipino et al.l 12009J) does 
not allow significant modifications of the underlying stellar 
yields. 

The SNIa rate for a simple stellar population (SSP) 
formed at a given time is cal culated assuming th e single 
degenerate scenario and the iMatteucci fe Recchil (|200ll ) 
Delay Time Distributio n (DTD). Th e convolution of this 
DTD with the SFR fsee lGreggicTl2005h gives the total SNIa 
rate. 

Stars - and baryonic processes at the galactic scale 
that need finer detail - are evolved between time-steps using 
sub-stepping of at least 1 Myr. During each sub-step, stars 
release mass and energy into the interstellar medium. In 
GallCS, the enriched material released in the late stages 
of stellar evolution is mixed to the cold phase, while the 
energy released from supernovae is used to re-heat the cold 
gas and return it to the hot phase in the halo. The rate of 
mass loss in the supernova-driven wind that flows out of 
the disc is directly proportional to the supernova rate (see 
below). 

Throughout this work, we assume chemical homogene- 
ity (instantaneous mixing), such that outflows caused 
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by feedback processes are assumed to have the same 
metallicity as the inter-stellar medium, though in reality 
the situation cannot be cap tured by our simple recipe 
l|Strickland fc Heckman| [2009) and newly produced metals 
are more likely to be ejected than the gas (see, for example, 
iRecchi et al.ll2004h . Note that, due to the fine sub-stepping 
used for the stellar evolution, ejecta from SNII and the 
contributions of single low-mass stars is implemented 
without the need to assume the instantaneous recycling 
approximation. 

Below we focus on the main free parameters of the model 
that we will change in order to study the build-up of the 
mass-metallicity relation. A complete, detailed description 
of Gall CS and a compendiu m of its entire features may be 
found in lHatton et al.l 



pendmi 
](2003) 



2.1 Initial Mass Function 

We use the Salpeter (|Salpeterlll955l l IMF 
4>(m) <x m -2 ' 35 



(1) 



in our simulations. The effect of changing the IMF is 
then investigated by replacing it with the Kennicutt 
l|Kennicutt et al.lll994h IMF, 



-1.4 
-2.5 ' 



m < 1M Q 
otherwise 



(2) 



in section [4] We take as the mass range 0.1 < m < 40M© 
for each IMF which determin es the normalisation. This is 
lower than adopted (65M@) bv lMaiolino et all (|2008l ); whose 
value would result in a higher O abundance on average of 
less than 0.1 dex. We note here that an even higher cutoff 
(greater than 80Mq) would lead to a change in the predicted 
O abundance of less than 0.3 dex, although we must rely on 
the extrapolation of the nucleosynthetic yields available in 
the literature, whereas the mass range m < 40Mq is where 
the O yields are more robust. We have performed tests that 
show, as expected, that changes in the yields effect only the 
normalization of the Mass-metallicity relation and not the 
slope. We therefore prefer to keep the same configuration as 
iPipino et all {2009), where the chemical evolution scheme 
(with a 40Mq upper limit) has been calibrated on obser- 
vations of the Milky Way. M oreover, there are indications 
(|Cescutti fc Chiappinil I2010T I that O production decreases 
with metallicity due to metal-dependent mass loss. In light 
of the above caveats, the reader should keep in mind that 
some fine-tuning of the normalization in the predicted mass- 
metallicity relation could be performed by acting on either 
the IMF or the O yields. 

2.2 Star Formation 

The star formation rate is given by 



SFR = q 



tdy 



(l + z 



,/3 



(3) 



where M co id is the mass of cold gas. The dynamical timescale 
tdyn is defined as the time taken for matter at the half- 
mass radius to reach either the opposite side of the galaxy 
(disc) or the centre of the galaxy (bulge) and is given in 



lHatton et ail ((2003). The free parameter a sets the efficiency 
of star formation and j3 controls its redshift (time) depen- 
dence (if any) . Th e GallCS fiducial v a lue of a (|Hatton et al.l 
20031) is 0.02 after iGuiderdoni et all l| 19981 ) and others (see 



Somerville fc Primack ( 19991 ) tables 4 and 5) have used val- 



ues in the range 0.01 < a < 0.25 with GallCS r estricted 
to the range 0.01 < a < 0.1 l|Hatton et al.ll2003l ). In this 
work we vary a over the range 0.01 < a < 0.05. There is 
observational evidence that the SFR was higher in the past 
llLillv et alJI 199rj; ISpergel et al.lll997l;lHelmboldt et al.ll2004l ; 
ljuneau et al. 20051 ; iFeulner et alj 20051 ) and so we focus on 
positive values of /3. Currently there is no ubiquitously ac- 
cepted theoretical explanation for this. One possibility is 
that the cosmic expansion means that galaxies were closer 
together in the past and could thus interact more easily than 
at later times giving a higher merger rate than at present. 
There is, however, evidence that the number of star-forming 
galaxies at z m 2 is much higher than the predicted num- 
ber of mergers l|Conrov et alj|2008t ) indicating that much of 
the s tar formation is not merger driven. Some (|Pipino et al.l 
2009) argue that positive feedback from the central super- 
massive black hole during its early growth period can ac- 
count for these observations and others that larger galaxies 
(that would have formed i nitially) have an over all higher 
cross section per unit mass jFerreras fc Silldl2003h and thus 
accrete gas faster than smaller galaxies formed more re- 
cently. Successfully simulating a predicted mass-metallicity 
relation that closely matches the observed one may eluci- 
date the mechanism by which galaxies have a higher SFR 
in the past. The GallCS fiducial value of f3 is 0, however, 
using f3 — 0.6 in orde r to incorporate the eff ect of rapid ac- 
cretion by cold flows, ICattaneo et al.l |2006) have managed 
to improve the fitting of the GallCS predicted luminosity 
function of Lyman-break galaxies to observations. Hence, 
we vary ft around this value. In section[5]we will discuss our 
predicted star formation rates in comparison with those that 
are observed in the galaxy samples with which we compare 
our predicted mass-metallicity relations. 



2.3 Supernova Feedback 

Massive stars will become type-II supernovae that release 
energy into the ISM causing a fraction of the gas to be 
ejected. If the energy is sufficient this ejected gas may leave 
the galaxy inhibiting star formation. This process is known 
as supernova feedback. W e model the outflow rate using the 
formula fromlsilkl l|2003l ) 



2£ ^SNI?SNSFR. 



(4) 



where SFR is given in equation [3] Here ??sn is the number 
supernovae per unit star forming mass, which we take as 
7.4xlO _3 /M and -Esn is the energy released by a single 
supernova assumed to have a constant value of 10 44 J. The 
escape velocity w e sc is given separately for bulges, discs and 
haloes in GallCS I (|Hatton et alj|2003h . The free parameter 
e determines the efficiency at which gas from supernovae is 
injected into the ISM. The mass-loading factor is then equal 
to the factors pre- multiplying the SFR and hence e _1 can 
be though of as the efficiency of mass- loading. The GallCS 
fiducial value is 0.3 and others have used values in the range 
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0.05 < g < 0.2 (see ISomerville fc Primackl l|l999l) tables 4 
and 5). ICattaneo et aL I (|2006l ) have found that e = 0.2 is the 
best-fitting value if GallCS is used to predict the luminosity 
function of Lyman-break galaxies. Hence, we initially vary 
e in the range < e < 0.3. 



2.4 Ejection of Matter from the Halo 

Type-II supernova and other galaxy forming processes may 
lead to the ejection of gas and metals from th e halo. These 
indiv idual processes are difficult to model l|Hatton et al.l 
2003) and the GallCS model accommodates them by stor- 
ing all ejected gas and metals ejected from the halo in a 
'reservoir' from which they may or may not be re-accreted 
at a later time. When matter is accreted from extra galac- 
tic sources a certain fraction is primordial and the rest is 
drawn from the reservoir. The efficiency of re-accretion from 
the reservoir is characterized by the parameter the halo 
recycling efficiency. If £ = then material ejected from the 
halo is lost permanently. If f = 1 then all matter accreted 
is drawn from the reservoir until it is depleted at which 
point any further gas accr e ted is assumed primordial. Re- 
cently, lOppenheimer et all (|2009l ) have studied the effects 
of recycling on the SFRs and stellar mass function of galax- 
ies in cosmological hydrodynamical simulations and have 
found that it is the dominant factor in galaxy growth at 
z < lj ejecta being the important factor at z > 2. Hence, in 
this work we focus on the supernova feedback efficiency and 
do not vary holding it constant at the GallCS fiducial 
value, 0.3, so that during any gas re-accretion process, 30% 
is drawn from the reservoir (provided it is not depleted) , the 
remainder being of primordial metallicity. 



3 DATA 

The gas metallicity mu st be determined usin g strong line 
metallicity diagnostics (|Maiolino et al.l I2008T ) which relies 
on the fact that the ratio of various strong optical emission 
lines depend upon the gas metallicity in a known manner. 
Thus these ratios must be calibrated against metallicity. 
Calibrations have only been performed in narrow metallicity 
ranges. These calibrations are often inconsistent with each 
other and can lead to different metallicity offsets of up to 
0.7 de x and the differenc e in th e shape of the curve is often 
large (|Kewlev fc EllisonI I2008T ). Since data measured at 
different redshifts may be measured using different optical 
lines it is essential to ensure a correct intercalibration 
between the data fits so that the correct evolution of the re- 
lat ion can be se e n. Th ese inter calibration issue s are t ackled 
by iNagao et~aT] (|2006l ) and by IMaiolino et al l (120081 ), who 
take a large sample of local galaxies spanning a wide range 
of metallicities (7.2 < 12 + \og{0/H) < 9.2, accurately 
measured by using both the electron temperature method 
and photoinoization models) and cross-calibrate the various 
strong line ratio diagnostics on the same metallicity scale. 



logAf* Kg 



0.07 11.18 
2.27 12.38 
3.54 12.87 



9.04 
8.99 
8.90 



Table 1 . Calibration constants used in equation \E\ as a function 
of redshift from IMaiolino et al.1 l|2008l l. 



where log Mo and Kq are free parameters that must be de- 
termined at each redshift to obtain the best-fitting to the 
observed data and are shown in Table [T] The calibration 
constants at z = 0.07 were derived using the data from 
I Kewlev fc EllisonI |2008T ). t he constan t s at z = 2.27 were 
found using the data from lErb et all (|2006T ) and the con- 
stant s corresponding to z = 3.54 were calculated using data 
from IMaiolino et"alT (|2008f ) . 



4 RESULTS 

In this section we test the effects of changing the parameters, 
in order to identify those that may improve the agreement 
with observations. A more quantitative discussion will in- 
volve only these parameters. Initially all qualitative results 
were obtained using the Salpeter IMF (equation [1} however 
the effect of changing to the Kenicutt IMF (equation [2} is 
investigated in section [ 



We t ake as our observed trend the AMAZE (M aiolino et al.l 
2008J) mass-metallicity relation 



4.1 Variation of Parameters 

4-1.1 Star Formation Efficiency 

Holding /3, e at their fiducial values (|Hatton et al.l [20031 ') . 
the parameter a was varied over the range 0.01 < a < 0.05 
to investigate the effects of changing the star formation 
efficiency on the predicted mass metallicity relation. Fig. 
[T] shows the output for galaxies at redshift 2.27 plotted as 
12 + log {O/H) Vs. log(M*/M @ ) where Af* is the stellar 
mass. Als o plot ted are the observations at z = 2.27 from 
lErb et all (|2006T ) as well as the calibrated best-fitting trend 
(equation [5]) at redshift 2.27 and 0.07. The corresponding 
plot for galaxies at redshift 3.54 is shown in Fig. [5] We note 
that, although the minimum baryonic mass is 2xlO 9 M0, 
it is possible to have galaxies whose stellar mass is less 
than this in the sample provided that their baryonic mass 
exceeds this minimum value. The plots show that increasing 
a has the effect of spawning a similar number of galaxies 
(at the same redshift) that have on average a higher mass 
and metallicity. Increasing a increases the SFR in direct 
proportion thus in the same Hubble time we have more 
stars formed and a larger proportion at a later stage in their 
life and so the stellar mass and metallicity is increased. 
We note from Fig. [2] that th e observed mass-metallicity 
relation IMaiolino et al.l 120081 ) would be best fitted using 
a low star formation efficiency for low-mass galaxies and 
a high star formation efficie ncy for high-m a ss gal axies. 
This supports the findings of IMaiolino et al.l (|2008h who 
argue that low-mass galaxies are characterized by low star 
formation efficiencies inhibiting chemical evolution. We will 
return to this issue later in section [5] 



12 + log {O/H) = -0.0864 (log M* - log M ) + K Q , 



(5) 



From a formal point of view, if we quantify the agreement 
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between model predictions and data in terms of the x 2 , 
we have values that monotonically decrease with increasing 
a, simply because the normalization of the predicted 
mass-metallicity relation increases and, on average, more 
model galaxies are in better agreement with the data. This 
trend, however, has the effect of predicting too many small 
galaxies at z ~ 3 that exhibit the metallicity of a typical 
z ~ 2 galaxy of the same mass. This trend is already 
present at a = 0.06 without inducing any improvement of 
the slope of the predicted relation. When discussing the 
yields, we showed that we adopt quite a conservative value 
for the O production, therefore we believe that values for 
a > 0.05 should be avoided even if they lead to a low 
value x 2 ■ Also, the reader should note here that during 
the preparation of the manuscript, several papers have 
introduced a more robust way to explore the parameter 
space and statistically handle the comparison between 
model prediction s and observat ions (see, for example, 
iBower et all l20ld ; iLu et alJlioioh . In particular, iLu et all 
1 20101 ) argue that the procedure adopted here (namely 
varying only one parameter at a time) does not allow one 
to uniformly explore the parameter space and that the 
"best-fitting by eye" values do not always coincide with the 
point that has the maximum likelihood in the parameters 
multi-dim ensional phase space. On the other hand, the 
procedure ILu et all l|2010l ) advocate may lead to a formal 
best-fit parameter set that is either unphysical or difficult 
to explain from the theoretical point of view. Therefore 
some priors on the parameters have to be adopted. The 
present work aims at probing the sensible range for some of 
those. 



4-1.2 Redshift Dependence of the SFR 

With a, e and £ held constant at their fiducial values j3 
was varied over the range < (3 < 1.25. The GallCS 
predictions for galaxies are shown in Fig. [3] (z ~ 2.27) and 
Fig- Bl ( z ~ 3.54). Inspection of the plots show that /3 has 
a distinguishable effect on the distribution of the mass and 
metallicities of the galaxies. At any given redshift increas- 
ing ft preserves (approximately) the number of galaxies 
predicted but the distribution favours more galaxies with 
a higher stellar mass and metallicity whilst preserving the 
trend which is a similar shape to the AMAZE calibration 
curve (equation [SJ. This follows from the fact that at 
constant redshift increasing j3 results in the term (1 + z) 13 
in equation being larger and thus it acts in a similar 
manner to a. The same number of mergers are predicted 
so that the mass distribution is similar but not identical 
since the SFR does play a role in deter mining the mass of 
individual galaxies (|Weidner et al.|[2007h . 
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Figure 1. Response of the predicted mass-metallicity relation to 
a change in the star formation efficiency for galaxies at redshift 
2.27. The blue crosses represent an individual galaxy with stellar 
mass and metallicity pr edicted by the m odel and the black aster- 
isks show the data from lErb et al. The green (upper) line 
shows the AMAZE best-fitting line for z = 0.07 and the black 
(lower) line is the same curve fitted to the data at z = 2.27. 



As in the previous section a simple \ 2 analysis would 
return lower results for increasing /3, whereas we believe 
that the values for fJ should not be higher than 1. 



4-1.3 Supernova Feedback Efficiency 

With q and /3 held constant at their fiducial values e was 
varied over the range 0.1 < e < 1.0. Following the discussion 
in section 12.31 this range is investigated in order to test 



the model however we do not use excessively large values 
in order to artificially fit relations hereafter. The plots in 



Figs. 5(a) and 5(b) show the mean mass and metallicity 
predicted in each of the six mass bins used by lErb et all 
(2006) at z = 2.27 and z = 3.54 respectively. Fig. |5(a) 
shows that the effect of increasing epsilon at z = 2.27 is to 
lower the average metallicity in each bin whilst preserving 
the overall shape of the relation. At higher values of e, 
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Figure 2. Response of the predicted mass-metallicity relation to 
a change in the star formation efficiency for galaxies at redshift 
3.54. The blue crosses represent an individual galaxy with stellar 
mass and metallicity pr edicted by the model and the black aster- 
isks show the data from lMaiolino et alj i2008h . The green (upper) 
line shows the AMAZE best-fitting line for z = 0.07 and the black 
(lower) line is the same curve fitted to the data at z = 3.54. 
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(b) = 1.25 

Figure 3. Response of the predicted mass-metallicity relation to 
a change in the redshift dependence of the star formation rate for 
galaxies at redshift 2.27. The blue crosses represent an individual 
galaxy with stellar mass and metallicity p redicted by the model 
and the black asterisks show the data from Erb ct al. (2006). The 
green (upper) line shows the AMAZE best-fitting line for z = 0.07 
and the black (lower) line is the same curve fitted to the data at 
z = 2.27. 



more gas is ejected from the galaxy thereby reducing the 
average metallicity in each bin. Thus the effect of changing 
e is to alter the offset of the relation without altering the 
slope. This effec t is co mparable to the one predicted by 
iFinlator fc Pavel |200ct ) with the difference that we have a 
self consistent chemical evolution that includes the finite 
lifetimes of both SNIa and SNII. At z = 3.54, Fig. [5(5)] 
shows that changing the value of e has very little effect 



on the relation, especially in the low mass regime, and 
basically preserves the shape of the relation. 



iMannucci et al.l (|2009l ). whose sample of LBGs at z ~ 3 
(discussed in the introduction), have found that the effective 
yield (the amount of metals synthesised and retained within 
the ISM per unit stellar mass) decreases with increasing 
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Figure 4. Response of the predicted mass-metallicity relation to 
a change in the redshift dependence of the star formation rate for 
galaxies at redshift 3.54. The blue crosses represent an individual 
galaxy with stellar mass and metallicity predic ted by the mode l 
and the black asterisks show the data from lMaiolino et al. (2008). 
The green (upper) line shows the AMAZE best-fitting line for 
z = 0.07 and the black (lower) line is the same curve fitted to the 
data at z = 3.54. 



naturally, regardless of the morphology, if the SSFR is larger 
in more massive galaxies and that galactic outflows are not 
needed to explain the relation. As noted in the introduction, 
none of the models that focus solely on different outflow 
processes only have been able to fit the observed relation 
at z = 3.54 and, taken together, these plots imply that 
outflows are only important in determining the low redshift 
relation whereas at higher redshifts, another mechanism is 
responsible for generating this relation, which may explain 
this lack of predictive power at higher redshifts. Considering 
these recent results, the findings of this section imply that 
it is the SFR-mass relation that generates the observed 
mass-metallicity relation and determines the slope at high 
redshift with outflows being important only for the low 
redshift properties. 



4.2 Variation of the IMF 

Using the free parameters held at their fiducial values the 
IMF was changed to the Kennicutt IMF (equation [2J . Fig. 
[6] shows the predicted relation for galaxies at redshifts 2.27 
and 3.54. Comparing Fig.[6]to Figs. I H4l it is clear that both 
IMFs predict similar shaped relations at both redshifts. As 
a matter of fact, the Kennicutt IMF predicts lower values 
than both Salpeter and observation at z — 2.27. At z = 3.54 
both IMFs predict similar distributions that are both in the 
same region as observed. Only the Salpeter IMF produces a 
distribution that has many galaxies with similar masses and 
metallicities to observations at both redshifts, which indi- 
cated that it may have been able to reproduce the observed 
relation given the right values of the free parameters. For 
this reason only the Salpeter IMF was investigated quanti- 
tatively to see if it could reproduce the observed relation at 
both redshifts. 



5 DISCUSSION 

Examination of the plots in Figs. [T] to 2] reveals that, no 
matter the value of any parameter, there is a definite trend 
in the predicted distribution of stellar mass and metallicity. 
This is the known trend of increasing metallicity with 
stellar mass that follows a shape similar to the AMAZE 
calibration curve (equation [5} at both redshifts. The plots 
also show a sc atter that is comparab le to the scatter seen in 
the data from iTremonti et al.l l|2004l ) at z ~ independent 
of the value of any of the free parameters. Finally most of 
the galaxies have met allicities that lie below the data from 
ITremonti et al.l (|2004l ) (z w 0) and the AMAZE best-fitting 
curve for z = 0.07 indicating that most of the galaxies had 
a lower metallicity in the past than today. This is con sistent 
with stellar and galactic chemical evolution theory (|Pagell 
1998) and is an indication that the model simulates the 
physics correctly. 



stellar mass for galaxies in their sample suggesting that 
galactic outflows cannot account for the shape of the 
mass-metallicity relation since their power is diminished in 
more massive galaxies and thus they cannot be responsible 
for the decreasing effective yields. Using chemical evolu- 
tion models f or ga laxies of varying morphological types 
ICalura et al.l ((2009) have found that the relation arises 



The results of section 14.11 showed that the model could 
reproduce the observed mass-metallicity relation and 
suggested that a good fit to the current observations could 
be achieved by a suitable choice of parameters and IMF. 
Section 14.1.31 revealed that the effects of varying e were 
much smaller than varying a and /3 hence only these two 
parameters were varied. Predictions were generated using 
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Figure 5. The change in the average mass-metallicity relation 
when e is varied at redshift 2.27 and 3.54. The purple circles, dark 
blue upward-facing triangles, light blue downward-facing trian- 
gles, green diamonds and yellow squares show the aver age mass 
and m etallicity in each of the six mass bins used by lErb et all 
(2006) when e = 0.1 , 0.3 , 0.5 , 0.8 and 1.0 respectively. The hori- 
zontal bars show the range of the mass bins and the vertical bars 
show one standard deviation in t he metallicity. T he black aster- 
isks and line show the data from lErb et all j2006t) (Fig. |5(a)| or 
iMaioIino et al.l l|2008h (Fig. |5(b)~[ | and the black (lower) line shows 
the AMAZE best-fitting line for the observed relations at the 
respective redshifts. The green (upper) line shows the AMAZE 
best-fitting line for the observed relation at z = 0.07. 



Figure 6. Predicted mass-metallicity relation using the Kenni- 
cutt IMF with the free parameters held at their fiducial values. 
The blue crosses represent an individual galaxy with stellar mass 
and metallicity predicted by the model and the green (upper) line 
shows the AMAZE best -fitting line for z = 0.07. The black aster- 
isks show the data fro m lErb et al.l (2006r ) (6(a)) or the data from 
IMaioIino et al.l l l2008h (6(b)) and the black (lower) line shows the 
AMAZE best-fitting line for z = 2.27 (6(a)) or z = 3.54 (6(b)) . 



parameters in the range 0.01 < a < 0.05 and < /3 < 1.25 
with e and £ held constant at their fiducial values. 



It was found that the parameters that best match the 
observed relations were a = 0.02 and /3 = 1.0. The plots 
in Fig. [8] show the predicted relations at both redshifts 
using the best-fitting parameters found in section [5] Fig. 
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shows that the observed relation is not predicted at 
z — 2.27, instead the predicted distribution has a large 
scatter and shows a similar s haped tren d to bu t lies mainly 



below the observations of lErb et all J200rj). Fig. |8(b) 



shows that the observed relation i Maiolino et al.l 120 08) is 
predicted at z = 3.54 with a comparable scatter. At both 
redshifts the m ajority of predic t ed ga laxies lie below the 
observations of iTremonti et all <|2004h and the AMAZE 
z = 0.07 best-fitting relation (equation O. At both redshifts 
we predict a flatter relation than observed. This effect is 
particularly clear at z = 2.27 where none of the predicted 
massive galaxies are consistent with the observed ones even 
taking into account the errors on the latter. The situation 
is less extreme in the high redshift case, nevertheless it may 
be a symptom of some missing physics in the model. We will 
thoroughly discuss such problems in the rest of this section. 
To gain better insights on the model galaxies it is worth, 
however, first having a look at the predicted SFRs, since 
the predicted flattening in the mass-metallicity relation 
may come from issues with the manner in which GallCS 
predicts them. Both plots show the predicted galaxies 
stratified by SFR into the three bins SFR < lOM yr _1 , 
fO < SFR < SOMoyr" 1 and SFR > 3OM yr" 1 . To some 
extent, the predicted SFRs should only be compared with 
each other and any direct comparison with experimentally 
observed values should be treated equivocally. In fact, 
GallCS predictions are average values over the entire 
galaxy and over the timestep, whereas observations look 
at more central regions (see below) and yield an instanta- 
neous SFR value. We can thus gain more insight from the 
comparison of the relative differences and evolution in the 
predicted and observed SFRs than their numerical values. 
Indeed, our predicted SFRs (s ee Fig. [7] below) are quite 
lower than those observed by IMaiolino et al.l (2008) who 
find an average SFR of ^lOOM ^yr -1 but are consistent 
with those from lErb" et al.l (|2006T ) who find a median value 
of ~ 3OM yr~ , although we do not reproduce some of 
their extreme cases (with SFRs ~ 3OOM yr _1 ). 

While a (small) fraction of the disagreement can be 
explained by a small difference in the adopted upper 
mass limit for the IMF between our model and the value 
adopted by the observational works, we believe that this 
is a gen eral problem of semi- analytic models (see, for 
example, Kh ochfar fc Silk] 120091 . and references therein). 
For instance, if the gas accretion and cooling are uneven 
over a timestep, the SFR may present spikes that are 
in better agreement wit h the values rep orted by both 
IMaiolino etail < |2008f ) and lErb et all (120061 ). The dynamics 
of these high redshift galaxies are not always consistent 
with a smooth star formin g disc, instead the y are highly 
turbulent (see, for example. iGenzel et al.ll2~008h and a large 
fra ction do not seem to rotate (for example, those observed 
by lLaw et all I2009T I . A similar disc ussion of the dynamic s 
of the AMAZE galaxies is given in iGnerucci et all l|2010l) . 
This star forming mode is not yet accounted for in GallCS 
and, perhaps, the high value for /3 that we obtain should 
be interpreted as a warning: the treatment of the star 
formation must be improved, possibly taking into account 
gas supplies that are enhanced by cold streams. Namely 
the required boosting in the SFR is not only given by an 
increase in the efficiency, but also by augmenting the gas 



supply. This has to preferentially happen in more massive 
galaxies, where the disagreement with observations occur 
and, as can be inferred from the results presented in the 
previous sections, a much better fit could be obtained by 
differentially increasing a and /3 at the high mass end with 
respect to their fiducial values, namely tho se that correctly 
fit the low mass end. On the other hand iMannucci et al.l 
(2009) have observed the SFR as a function of log (M*/M ) 
for another sample of galaxies at z ~ 3 and have found, 
on average, a lower SFR than in the AMAZE sample, 
while at the same time the O abundances are similar 
to those we are using in this work. In particular, in the 
mass range 9.5 < log (M„/M ) < 10.5 the average SFR 
that we predict is (SFR) = 3O.OM yr _1 at z = 2.27. We 
note that this value has been corrected to account for the 
difference in IMFs. The range of the deviation from this 
average is (approximately) 75M yr~ at the three standard 
deviation level. At z = 3.54 and for similar masses we find 
(SFR) = 9.4 ± 26.4M yr _1 at the three standard deviation 
level. 

We note that the re are many m o re pr edicted galaxies 
in our bin than IMannucci et al.l |2009T ) and that the 
standard deviation in our SFR is due to this large number 
of predicted galaxies whereas theirs is due to experimental 
uncertainties involved with their observation. Therefore, an- 
other possibility, is that, while AMAZE probes the systems 
with the highest SFR, the model adopted here gives values 
that are more "typical" for high redshift galaxiefl. Never- 
theless, a successful model should predict also the extreme 
galaxies, unless those with extreme SFRs cannot be cap- 
tured by either the spatial resolution or the lack of physics. 
Finally, we note that the observed SFRs discussed above 
have been d erived by assuming e xponentially declining 
SF histories. iMaraston et al.l (|2010f l have shown that this 
assumption might lead to an overestimate of the SFR with 
respect to other, more realistic SF histories, especially when 
the age is not constrained to be larger then the e-folding 
time of the assumed e xponentially declining SFR. Similarly, 
IMaiolino et al.l (|2008T l have noted that the fits that yield 
the SFR in the AMAZE sample sometimes return unrealis- 
tically small ages (below 50 Myr) if not suitably constrained. 

In Fig. we plot the predicted SFR against stellar 

mass, ft is evident from this Fig. (as well as Fig. [S] discussed 

below) that galaxies with larger stellar mass and metallic- 

ities have higher SFRs and thus our predicted relation is 

one between SFR and mass, consistent with the recent pre- 

i 1 i — — i 

dictions of Calura et al. (2009) as well as the observational 

findings of IMannucci et al.l (|2009j). It is also clear from the 

plots that we do find that, at any given mass, more pristine 

galaxies have a large SFR and viceversa, in agreement with 

the recent suggestion by IMannucci et al.l (|2010f l that the 

mass-metallicity relation is, in fact, a mass-SFR-metallicity 

surface. Howeve r, if we use the new "variable" suggested by 

IMannucci et~al] (|2010l ). namely log (M») - 0.32 log (SFR) 

instead of log(M*), as the abscissa then we do not find a 



1 The reader should note that, given all these issues, we chose not 
to apply any cutoff based on the SFR to our prediction. Instead 
we show all the star forming galaxies in each bin. 
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Figure 7. The predicted SFR as a function of stellar mass for 
galaxies at redshift z = 2.27 and z = 3.54. The red asterisks show 
individual predicted galaxies with a given stellar mass and SFR. 



significant decrease in the scatter. Here we do not make any 
infer ences or draw any conc lusions about the relationship 
that iMannucci et al.l (|2010l ) have proposed but leave a 
more comprehensive comparison for future work where we 
intend to use a forthcoming, updated and refined version 
of GallCS that is better calibrated on the z = 2.27 and 
2 = 3.54 mass-metallicity relation. 

Our inability to fit the relation at z = 2.27 may be due 
to differences between the regions of the galaxies sampled 
by observ ation and by Ga lfCS. For instance, the obser- 
vations of lErb et all (20061 ) are obtained using a long slit 



spectrometer with an aperture of 0.75", corresponding to 
a radius (half-aperture) on th e galaxy of only 3.6 kpc at 
z = 2 27 ijMaiolino et al.ll20Qgl ). This implies that lErb et~ai1 
(2006) only sample the central region of the most massive 
- hence larger - galaxies, whereas less massive systems are 
probed out to larger radii. Local galaxies have metallicity 
gradients that, combined with the fixed aperture, bias the 
mass-metallicity relation. As shown by Cresci et al. (20f0, 
in press) high-redshift star forming systems have a complex 
metallicity structure, however further study is needed to as- 
sess its impact on the observed relation. On the theoretical 
side, GallC S calculates the av erage metallicity over the en- 
tire galaxy l|Pipino et al.l [20091 ). including the outer regions 
which are metal poor. We will come back on this aspect of 
this problem later when discussin g the role of the bulge and 
the burst components. AMAZE (|Maiolino et all 120081 ) use 
integral field spectroscopy to measure their metallicities, 
extracting the spectral information using a 0.7 " aperture 
which corresponds to a similar radius to the Erb et al.l 
(2006) sampl e. Moreover, the observed size evolution (see, 
for example, lLaw et all 2007) of LBGs between z ~ 3.5 
and z ~ 2, is r ather mild, therefore we deem it unlikely 
that SEE et al.l |2006l ) sample a much smaller portion of 
these galaxies - at a given fixed mass - than AMAZE, in 
order to justify an even larger departure of the models 
from the observations at z — 2.27 than at z — 3.54. Whilst 
these effects may be present (indeed we cannot rule them 
out) and can partly explain the difference between our 
models' predictions and observations, we believe that the 
disagreement at the high mass end has more to do with the 
treatment of the physics and the assumptions in the model. 

To better explain the lack of chemical evolution in the most 
massive galaxies we have performed several tests u sing a 
chemical evolution model (|Pipino fc Matteuccil 120041 ) with 
the same setup (in terms of both IMF and stellar yields) and 
allowing for different infall and star formation timescales. 
These models, that follow the evolution of a single galaxy, 
have the advantage of giving a self-consistent description of 
the chemical evolution, whilst using a simpler parametrisa- 
tion for the accretion and the star formation histories than 
we use in GallCS. In passing, we note that these models 
are, in essence, the same ones that are used bv lCalura et al.l 
(2009) to study the AMAZE galaxies in the context of their 
monolithic collapse. As far as the adopted star formation 
rate is concerned, it is given by the same formula as GallCS 
(equation [3]), with the absence of any redshift dependence 
i.e. P = 0. The star formation timescale isF = tdyn/ct and 
we assume an exponentially decreasing gas infall law (i.e. 
oc e _t,/r , where r is the infall timescale), that is not linked 
to any Dark Matter accretion history. By running these 
chemical evolution models for a large range of different 
combinations of r and igp, we confirm that a relatively high 
O abundance (12 + log (O / H) > 7.5) is attained within a 
few 100 Myr. At the same the time stellar mass is quickly 
built up. This is why GallCS, whose predicted galaxies 
feature a differing SFR and different accretion/merging 
histories, predicts a steep mass-metallicity relation at 
low masses/younger (relatively) ages. What differs is the 
subsequent evolution. For a star formation timescale much 
shorter than the infall timescale, the chemical evolution 
model show a metallicity 12 + log (O/H) > 8.5 already at 
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200 Myr and a very mild evolution after that. On the other 
hand, when the star formation timescale is comparable to 
the infall timescale and both are of the order of 1 Gyr, 
galaxies attain a metallicity 12 + log (O/H) « 8 in the same 
timescale and a further 0.5 dex increase in the O abundance 
would take more than 1 Gyr (i.e. the time elapsed between 
z=3 and z=2). In other words, since this latter case is 
typical in GallCS, it is easy for the model galaxies to 
quickly attain 12 + log (O/H) ~ 7.5 — 8 at redshift 3, 
hence reproducing the average observed value, but then 
the predicted metal enrichment in the redshift range 3-2 
is slower than that inferred from the observations because 
the star formation is not efficient enough. The former case 
(star formation timescale shorter than infall timescale), 
ins tead, charact e rises t he monolithic collapse models used 
by ICalura et all |2009l ) to reproduce the AMAZE data. It 
must also be remembered that the /3-boosting at z=3 is by 
a factor of 4, whereas the boosting at z=2 is by a factor 
of 3. Hence, the redshift dependence of the SFR in such a 
"narrow" range is quite mild and will not necessarily drive 
major changes in the predicted mass and metallicities (see 
Section 4) in this time interval. 

Moreover, at higher masses, a fair fraction (approximately 
20%) of model galaxies have had at least one merger. 
Therefore some older galaxies could have become quiescent, 
or at least have a sizable number of stars in the bulge 
component. In this latter case, even if the O abundances 
are quite high (12 + log (O/H) > 8.5) in the gas associated 
with the bulge component, the SFR and the gas mass in 
this component are much lower than in the disc. Were 
such a strict separation between disc and bulge in the 
models present also in actual galaxies, the bulge metallicity 
would not be easily observed. This seems to be at variance 
with actual observations (see, for example, lLaw et al.l 
12007ft that do not detect a drop in the rest-frame UV 
flux in the central regions of the galaxies. This fact may 
hint that the disc/burst /bulge decomposition of model 
galaxies that well characterises local galaxies is not a good 
approximation of the complex structure of high redshift star 
forming systems. Modifications such that, for instance, the 
metal rich gas from the bulge component could enrich the 
discs, may help in reconciling such models with observations. 

At the same time, as it can clearly be perceived from 
the fact that there are many more galaxies in the z=2 
figures as opposed to the z=3 ones, we predict a large 
number of new galaxies that appear during the elapsed 
time interval. The youngest ones have stellar masses below 
10 10 Mo and low metallicities, thus biasing the predicted 
mass-metallicity relation to lower values and affecting the 
predicted evolution. Therefore, in order to improve the 
agreement with observations, one can apply a SFR/age cut 
to the predicted sample of galaxies. Given the difficulties in 
predicting the correct SFRs - see above - this has not been 
attempted. 

A closer inspection of our galaxies show that the metal- 
licity increases with decreasing gas mass fraction almost 
as expected in a simple closed box model of chemical 
evolution, namely Z oc ln(M t ot/M gas ). The gas fraction 
is almost constant (around 0.8) with total baryonic mass 



(at a fixed redshift) in the high mass regime, whereas it 
exhibits a large scatter at the low-mass end. On the other 
hand it strongly decreases with stellar mass, although with 
increasing scatter. This is because the SFR is constant for 
galaxies with stellar mass below 1O 1O M0, becoming pro- 
portional to the stella r mass above this limit (c.f. Fig. 9 in 
iMannucci et al.l (|2009l )V Only for the most massive galaxies 
(in terms of their stellar mass), does the gas fraction rise 
again to (approximately) 0.8. This is basically what we 
see in the mass-metallicity plots: the metallicity increases 
with stellar mass up to ~ 1O 1O M0 where it then flattens 
out in the most massive galaxies (where the gas fraction 
is higher again). As explained above, in these systems the 
star formation is not efficient enough com pared to the 
inflow of primordial gas. Indeed, as shown by ICalura et al.l 
(2009), the high mass end of the z=3 mass-metallicity 
relation cannot be reproduced by the "local" disc, but 
proto-spheroids with very high SFR, such that their O 
abundance quickly jumps to solar values in a few Myr, 
are needed. At variance with such a model, where the 
galaxy morphology is assumed, GallCS galaxies may have 
three components ( disc, bulge and burst, see Section 2 and 
lHatton et all [2003), for more details) that co-exist at the 
same time and whose presence is linked to the evolutionary 
path of the galaxies. Here we remind the reader that we 
show the abundances predicted for the disc component of 
each galaxy only in order to have a fair comparison with the 
observations that sample quite a large region compared to 
the assumed sizes of both the bulge and burst components 
in the model (whenever they are present). The results do 
not change if we show the abundances averaged over the 
three components because the mass in the cold phase of the 
disc is much larger than that in the bulge (where we predict 
gas mass fraction well below 0.1) or in the burst. In these 
latter components we predict 12 + log (O/H) > 8.5 in the 
majority of the galaxies, however, such a high O abundance 
is diluted in the mass averaging. We have already noted 
that observations sample the inner regions of the galaxies. 
On the other hand, we present the results pertaining to 
the entire disc because the nominal spatial resolution of 
GallCS is ~30 kpc. If we had enough spatial resolution to 
make predictions about the inner ~4 kpc only, it is likely 
that: i) the local density would make the local SFR higher, 
hence the metal enrichment quicker and ii) the dilution of 
the metals in the bulge and in the burst by means of the 
inner disc gas would be lower (meaning a higher metallicity 
in that region of the disc, lower fractional contribution in 
the mass averaging). Finally, we note that on the abscissa 
of all the Figs, discussed so far, we plot the total stellar 
mass, summed over the three component. Since bulges and 
bursts tend to be more frequent at the high mass end, 
where galaxies are older and had more time to evolve and 
merge, this implies that we stretch the mass axis without an 
increase in the O abundance, hence obtaining an artificially 
flatter (by a small amount) and more extended to high 
masses mass-metallicity relation than the one expected for 
a pure disc. 

One possible way out would involve making predic- 
tions that take into account the three component, without 
computing an average O abundance. For instance, one 
could use the discs to explain the low mass end of the 
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mass metallicity relation, whereas bursts and bulges could 
explain the high mass end (see also lCalura et al"1l2009h . As 
noted above, the bulges cannot be directly compared with 
the observations because it is assumed that the SFR in 
this component can involve only gas returned from stars, 
hence the predicted SFR are lower than those in the disc. 
Using the bulges instead of the discs at the high mass end, 
where they are more abundant, will lead to an increase in 
the slope of the predicted mass-metallicity relation, at the 
expense of an even poorer agreement between the observed 
and predicted SFR. 

At z—2, nearly 15% of the galaxies feature a burst 
component. In such a component, the galaxies exhibit SFRs 
comparable to, or even a factor of ~2 higher than, the disc 
and the predicted O abundances are 8 < 12+log (O/H) < 9, 
so one may be tempted to use only the predicted properties 
of the burst only to compare with observations. The 
problems with this last scenario are numerous. In the 
first place, the burst component is assumed to have a 
radius that is typically below 1 kpc, hence smaller than 
both that of the disc and the aperture used by observers. 
We find it difficult to explain the properties of observed 
galaxies only by means of such a centrally concentrated 
component. Moreover, one has to devise a reason to bias 
the observations in favour of a "burst" only above a given 
mass (~ 1O 1O M0) in order to steepen the mass-metallicity 
relation at the high mass end. In part, such a bias can 
be granted by the fact that observed samples are selected 
according to their SFE0. However, since the model switches 
the burst component on as a consequence of mergers, the 
hypothetical solution of having a greater frequency of burst 
components in the more massive galaxies is at variance 
with the fact that a fraction of the observed galaxies do not 
show such signatu res. Moreover, recent observa tions (for 
example those of iForster Schreiber et al.l I2009T ) indicate 
that galaxies that are observed with a regularly rotating 
disc are preferentially higher mass systems. However, these 
discs still have large amounts of random motions and 
turbulence can be fairly thick. Future developments of 
the models should implement the possibility that "bursts" 
are not only centrally located. Instead SF clumps would 
be distributed in the disc, as also discussed above when 
comparing predicted and observed SFRs. 

The net effect of the changes suggested above should create 
a typical galaxy formation path in which, broadly speaking, 
a increases with galactic mass. Indeed a direct increase of 
a with mass m ay be achieved by directly incorporating the 
suggestion by |Pipino et all 12009) on the positive feedback 
by the central black hole. Such a change can be easily 
implemented in monolithic formation models, however it is 
more difficult to devise a mechanism in a scenario where 
larger galaxies are formed through mergers of smaller 
companion. This is where the study of the mass-metallicity 
relation at z ~ 3 provides leverage: many (predicted) 
galaxies have not undergone any mergers yet and sot he 
effect of these mergers is kept at a minimum level. These 



2 At least at z=2. In fact, the predicted SFRs in the burst com- 
ponents at z=3 are still lower than the observed ones. 



mergers tend to produce high mass galaxies with lower 
metallicities than those that have evolved in the standard 
manner and thus (in general) flatten the relation. Hence 
the observed mass-metallicity relation is, in fact, giving us 
insight into the SFR-mass and inflow-mass relations. To be 
more quantitative, from the investigation done in Section 
4 we can infer that if the value of a is kept at the fiducial 
GallCS value at the low mass end, whereas its value is 
increased by a factor of at least 3 at the high mass end, the 
predicted mass-metallicity relation would be considerably 
steeper and in better agreement with observations. The 
reader should note that this is an illustrative example to 
provide insight into the manner in which current models 
should be modified in order to improve them in this respect 
and that unnaturally fixing a as a function of stellar mass is 
only an artificial solution and is not the correct method by 
which this can be achieved. In passing, we note that these 
solutions may also alleviate the problems that semi-analytic 
models have in reproducing the abund ance ratio-mass re - 
lations in present-day elliptical galaxies l|Pipino et alj|2009h . 

We also note that, although variation of the parameters 
significantly changes the prediction of the mass-metallicity 
relation at high redshift, the effect at low redshift is much 
less evident. 

On the other hand, examination of the predicted distribu- 
tion at z = shows that the evolution in metallicity between 
z = 3.54 and z = 0.07 agrees with the ob served evolution 
jTremonti et alj|2004 IMaiolino et alJl200S( ) . This is shown 
in Fig . [51 Previous models ( for example |Pe Lucia et al.l 
l|2004 l. iKobavashi et all \200H) and others discussed in the 
introduction) often fit the relation at one redshift only 
and the redsh i ft evo lution is not tested to z = (see 
IMaiolino et al.l (|2008h section 7.5 for a discussion). We have 
fitted the observed relations at both z = 3.54 and z « 
and hence the model is more tightly constrained. We have 
not been entirely successful in fitting the observed relation 
at all redshifts however we have attempted to constrain it 
in different regimes and predict the observed evolution of 
the relation from z = 3.54 to z = 0. Furthermore, we have 
found that outflows are not the origin of the relation and 
only effect its low redshift evolution and hence we have not 
needed to impose different outflow models to fit at specific 
redshifts. Nor have we needed to fine-tune the efficiency of 
supernova feedback. 



6 CONCLUSIONS 

The relation between stellar mass and gas-phase metallicity 
has been well documented over the range < z < 3.5 4 
jTremonti et alj2004 lErb et alfeOOrj ; IMaiolino etalll2008l ) . 
Despite this many hydrodynamical and semi-analytic N- 
body simulations of galaxy formation have been unable to 
predict the correct relation over a suitab l e redshift range 
ilDe Lucia et alj|2004 iDe Rossi et~ai1l2007l ; iMouchine et all 
2008). In this paper we have fitted the observed relation 
at z ~ 3.54 u s ing a n updated version of GallCS given 
in iPipino et al.l (|2009T i and have reproduced the observed 
relation at z ~ 0. Th e model uses a AC DM cosmology 
taken from WMAP3 l|Spergel et all 120071 ') and includes 
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Figure 8. Predicted mass-metallicity relation using the best- 
fitting parameters a = 0.02, = 1.0, e = 0.15 and £ = 0.3 
and the Salpeter IMF (equation [2j . The yellow crosses represent 
galaxies with SFR < lOAfQyr -1 , the light blue triangles repre- 
sent galaxies with 10 < SFR < 3OJV/0yr~ 1 and the dark blue 
circles represent galaxies with SFR > 3OM0yr — . The green (up- 
per) line shows the AMAZE best-fittin g line for z = .07 and 
the black asterisks show the data from [Erb et alj {2006 ) | |8(a) \ 
or lMaiolino et al.l feOOSl ) j8(b)| l. The black (lower) line shows the 
AMAZE best-fitting line for z = 2.27 (8(a)) or z = 3.54 \8(b)\ . 
The red crosses show the average m a ss an d metallicity in each 
of the 6 mass bins used by Erb et al ] l|2006l l. the horizontal bars 
show the mass bins and the vertical bars show one standard de- 
viation in the metallicity. 
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Figure 9. The predicted evolution from z = 3.54 to z 0. The 
black crosses show the a verage mass and met allicity in each of 
the six mass bins used bv lMaiolino et al. I J2008h at z = 3.54 using 
the best-fitting parameters. The horizontal bars show the range of 
each mass bin and the vertical bars show one standard deviation 
in the metallicity. The AMAZE best- fitting curve (equation [5J at 
z = 0.07 is shown by the green line and the blue crosses show the 
predicted mass-metallicity relation at z fa 0. 



outflows due to type -II supernovae since previous work 
(|De Rossi et alJ 120071 ) has shown that these are needed 
to reproduce the observed relation. Many free parameters 
are included that control the star formation efficiency, 
redshift dependence of the SFR, type-II supernova feedback 
efficiency and the halo recycling efficiency which were 
varied separ ately in order to achie ve the best possible fit to 
observation (|Maiolino et al.l I2008T ) . Here we summarise the 
results of the simulations. 

Firstly, we have found that the effect of varying the 
supernova feedback efficiency has little impact on the rela- 
tion at high redshift but does have a small effect at lower 
values. We have also predicted that more massive galaxies 
have higher SFRs than less massive ones. These two results 
taken together support recent findings (Ma nnucci et al.l 
l2009l ; ICalura et alj|200gh that it is the SFR-mass relation in 
galaxies and not galactic outflows that is responsible for the 
origin and evolution of the relation although the cumulative 
effects of outflows do affect the low redshift evolution. 
Secondly, a better agreement between the predicted and the 
observed average metallicity at a given redshift can be ob- 
tained by assuming that the SFR that has a strong redshift 
dependence, proportional to 1 + z, which is slightly stronger 
than other models have used in the past (jCattaneo et al.l 
(2006) uses (1 + z) 0,6 ). At the same time, the predicted 
SFRs are lower than the observed values at z ~ 3, but are 
consistent at lower redshifts. These facts may point to the 
need for a revision of the SFR recipes in a future generation 
of semi-analytical models. Thirdly, the observed relation at 
z — 3.54 is reproduced well with a scatter in the distribution 



The Mass-Metallicity Relation at High Redshift 15 



of stellar mass and m etallicity comparable to observations 
ijMaiolino et al.l l2008f ). However, the predicted relation is 
flatter at the high m ass end. Also, w e fail to reproduce the 
relation at z — 2.27 (lErb et al.l 120061 ). where the flattening 
in the high mass regime becomes more evident. We discuss 
several reasons for this disagreement, stemming from both 
theoretical and observational biases. We argue that if 
observations are preferentially selecting galaxies with high 
SFRs and the measured abundances mirror those in these 
"bursts" rather than averages over the three components, 
we might solve the problem of the flattening by using the 
abundances predicted in the "burst" rather than in the 
disc component of the model galaxies. This is a rather ad 
hoc solution and presumably means that future models 
should take into account that, at high redshifts, either 
the SF occurs in clumps rather than in a smooth disc or 
the star formation efficiency is boosted preferentially in 
the mos t massive galaxies by, for example, positive AGN 
feedback Iplpino et all (|2009l ). We note, that the model does 
reproduce the observed evolution to z = 0. The efficiency 
of outflows was found to have only a minimal effect on the 
predicted distribution. Finally, we have seen that low-mass 
galaxies are best fitted to the relation using a low star 
formation efficiency whereas the converse is true regarding 
hi gh-mass galaxie s . This observation supports the findings 
of iMaiolino et all ([2008 ) and thus galaxy downsizing may 
be the origin of the ubiquitously observed relation. 
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